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Abstract. We examine the most recent observational con- 
straints arising from i) small-scale and large-scale Galactic dy- 
namical properties, ii) star counts of population I and II stars at 
faint magnitude and iii) microlensing experiments towards the 
Large Magellanic Cloud and the Galactic centre. From these 
constraints, we determine the halo and disk stellar mass func- 
tions and stellar content down to the bottom of the main se- 
quence, which yields the normalization of the halo/disk total 
stellar population, and we infer the contributions of sub-stellar 
objects to the mass budget of the various Galactic regions. 

The consistent analysis of star counts and of the overall 
microlensing observations in the Bulge are compatible with a 
small contribution of brown dwarfs to the Galactic mass bud- 
get pbd/p* < 0.2. However the separate bulge/disk analysis 
based on the bulge clump giants is compatible with a substan- 
tial population of disk brown dwarfs, £bd/£* < 1. More 
statistics of microlensing events towards the Galactic center 
and a better determination of the velocity dispersions in the 
bulge should break this degeneracy of solutions. 

For the halo, we show that a steep mass-function in the dark 
halo is excluded and that low-mass stars and brown dwarfs rep- 
resent a negligible fraction of the halo dark matter, and thus 
of the observed events towards the LMC. The nature of these 
events remains a puzzle and halo white dwarfs remain the least 
unlikely candidates. 

Key words : stars : low-mass, brown dwarfs — The Galaxy 
: stellar content — The Galaxy : halo — Cosmology : dark 
matter 



1. Introduction 

There is compelling evidence for believing that most of the 
matter in the Universe is under the form of dark, yet unob- 
served components. The observed density of baryons in galax- 
ies f^star+gas ~ 0.003 hr 1 (h is the Hubble constant in units 
of 100 km.s _1 Mpc _1 ) may represent only a fraction of the 
value predicted by primordial nucleosynthesis, CIb ~ 0.01 h~ 2 
(Copi, Schramm & Turner, 1995). There is also evidence that 
spiral galaxies are surrounded by a large amount of non- 
luminous mass of unknown nature. These two facts suggest 



that baryonic dark matter is a possible candidate for halo dark 
matter. A breakthrough in this longstanding, unsolved prob- 
lem has been accomplished recently with the developement of 
microlensing experiments, by inferring the presence of dark 
objects through their gravitational effect on luminous matter. 
The analysis of the EROS (Aubourg et al., 1993; Ansari et 
al., 1996) and MACHO (Alcock et al., 1993; 1996) observa- 
tions towards the LMC, in particular, reveal the presence of 
some baryonic matter in the Galactic halo, although the inferred 
optical depth shows that this dark baryonic matter probably 
does not provides all the sought missing mass. On the other 
hand, microlensing experiments towards the central regions of 
the Galaxy yield a mass density of unseen star-like objects 
about three times larger than the value expected from standard 
disk+bulge models (Udalski et al. 1994; Alcock et al., 1997). 
These two results from microlensing experiments - the lack of 
events in the halo and the excess of events in the bulge - lead 
to the tempting conclusion that more (resp. less) galactic dark 
matter than expected previously resides in the disk (resp. the 
halo). On the other hand, severe constraints on the amount of 
dark matter in the disk and in the halo arise from the observed 
large-scale kinematic properties, i.e. the rotation curve of the 
Galaxy at distances larger than the observed luminous distri- 
bution, as well as from the small-scale dynamical properties, 
i.e. the stellar velocity dispersion in the solar neighborhood. 
This latter information is derived from the measurement of the 
vertical acceleration due to the galactic potential near the Sun, 
which yields the determination of the local dynamic surface 
density. Finally, star counts have now been obtained at very 
faint magnitudes, either from the HST or from ground-based 
deep-magnitude surveys, and provide stringent constraints on 
different parameters entering galactic modeling such as scale 
lengths, scale heights, and more indirectly stellar mass func- 
tions. 

These two latter type of constraints, galactic dynamics and 
star counts, have usually suggested that most of the Galactic 
missing mass resides in the outer so-called dark halo. There 
is then an apparent conflict between the afore-mentioned mi- 
crolensing results and standard astronomical observations. This 
paradigm stresses the need to reconsider the standard model of 
the Galaxy, where all dark matter resides in the halo, and to 
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derive a Galactic mass-model consistent with the three types 
of observational constraints, kinematics, star counts and mi- 
crolensing. This is the aim of the present study. 

The calculations will be presented in two joint papers. 
In the present one (Paper I), we examine in detail all obser- 
vational constraints arising from the most recent determina- 
tions of small-scale and large-scale Galactic dynamic prop- 
erties (§2), from star counts in the Galactic disk, bulge and 
spheroid (§3), and from microlensing experiments towards both 
the Large Magellanic Cloud (LMC) and the Galactic center 
(§4). A consistent analysis of all these observations yields the 
determination of the stellar mass functions, slope, normaliza- 
tion and minimum mass, and thus of the amount of mass under 
the form of stellar and sub-stellar objects in the various regions 
of the Galaxy, disk, bulge, spheroid and dark halo. 

In paper II, these results will be used to derive a consistent 
model for the Galaxy, confronted to all types of observational 
constraints. We will discuss in detail the two possibilities of 
an essentially non-baryonic and a dominantly baryonic mass 
model for the Galaxy (Mera, Chabrier & Schaeffer, 1997). 

2. Kinematic constraints 

In this section we examine the most recent observational de- 
terminations of the kinematic properties of the Galaxy. We first 
consider the asymptotic circular velocity, in the Milky Way and 
in other spiral galaxies, and then the velocity distribution in our 
Galaxy up to 60 kpc. Different determinations of the total mass 
of the Galaxy, converted into equivalent circular rotation veloc- 
ities v 2 rot = GM(r)/r, are also examined. In the second part 
of this section, we focus on the determination of the local sur- 
face density, inferred from the measurement of the gravitational 
acceleration of the Galactic potential near the Galactic mid- 
plane. From all these observations, we derive the most likely 
values for the Galactic mass, rotation velocity, and amount of 
dark matter near the Sun, to be reproduced by the final galactic 
model. 

2.1. Large scale velocities 
2.1.1. Circular velocity 

The most recent determinations of the circular velocity as a 
function of the Galactic radius have been reviewed by Fich and 
Tremaine (1991). Their compilation of the measurements of the 
outer rotation curve from CO, HI and HII observations yields 
the accurate determination of a velocity of about 220 km.s -1 , 
constant within ^ 10% up to 14 kpc, with extremely small error 
bars (a few %). At larger distance there is a hint for a small in- 
crease of the rotation curve, but the uncertainty in the data does 
not allow any robust conclusion. In any case, there is no sign 
for a significant decrease of the velocity up to about 20 kpc. 
More recent observations of neutral hydrogen out to 2.5 Rq, 
where R Q ~ 8 kpc is the galactocentric position of the Sun, 
yield a slightly lower value t> rot = 200 ± 10 km.s -1 (Merri- 
field, 1992). On the other hand, Schechter et al (1989) report 



v mt = 248 ± 16 km.s -1 from the kinematics of carbon stars 
in the outer Galaxy. As noted by Kuijken & Tremaine (1994), 
a difference between the value determined from stellar and HI 
tracers is expected, because of the ellipticity of the galactic po- 
tential. This will decrease the afore-mentioned stellar value. 

The determination of the rotation curves of external galax- 
ies (Casertano and Van Albada 1990) characterizes what these 
authors call "bright galaxies", with circular velocities bracket- 
ted between 180 and 260 km.s -1 . These rotation curves show 
occasionally a slight drop beyond 20 kpc, which for all of them 
remains within less than 20% at 35 or 40 kpc, almost within 
the error bars. The rotation curves are thus flat within the error 
bar determination (<~ 10 — 20%) up to the furthest distance at 
which hydrogen is detectable, i.e more than twice the radius 
of the visible stellar component. Although no rotation velocity 
measurements exist for our Galaxy beyond 19 kpc, Fich and 
Tremaine (1991) argue that, by comparison with other galax- 
ies, our rotation curve should extend similarly up to at least 35 
kpc. 

Note that the Galactic rotation is well confined to the 
disk plane. At 1 kpc above the disk, metal-rich stars show 
similar rotation velocities and are considered as belonging to 
the (thick) disk population. Most of the metal-depleted stars 
([Fe/H] < —1.5), which probe essentially the halo popula- 
tion, rotate only with <~ 30 — 40 km.s -1 assuming a rotation 
velocity 220 km.s -1 for the Local Standard of Rest, arguing 
for a non-rotating, or slowly rotating halo (Beers and Sommer- 
Larsen, 1995). 

2.1.2. Velocity distribution up to 60 kpc 

Attempts to determine the gravitational potential of our Galaxy 
at large distances, from the study of the motion of globular 
clusters, require some modelling of the velocity distribution in 
the Galactic halo, which enters the Jeans equation (Binney & 
Tremaine 1987). The compilation of Harris and Racine (1979) 
led Frenk and White (1980) to use a 66 cluster sample, with 
various assumptions for the shape of the Galactic halo, to deter- 
mine the mass of the Galaxy within 33 kpc. These clusters show 
an overall 60 ± 26 km.s -1 circular rotation (assuming vi sr = 
220 km.s -1 ), again advocating for a slow rotation of the halo, 
but also for important non-circular motion. Translating the dy- 
namically measured mass of the halo within a given radius into 
an effective circular rotation velocity (v^ ot = GMj R), their 
results can be summarized as follows : this velocity always lies 
within the range 200 < v < 319 km.s -1 , for all models, at the 
90% confidence level. With the additional constraint that the 
globular cluster distribution is not more flattened in the outer 
region than in the inner region, the allowed velocity-range re- 
duces to 200-265 km.s -1 . A second restriction stems from the 
fact that if the velocity distribution is not isotropic, the radial 
rms velocity a r must be larger than the tangential velocity v t , 
because of infall towards the centre of the Galaxy. This leads to 
the final range 224km.s -1 < w rot (33kpc) < 265 km.s -1 , the 
uncertainty being due to systematic errors. 
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The more distant the constraint, the more relevant it is to de- 
termine the dark matter content of the Galaxy, but the more un- 
certain the data. The sample of 17 globular clusters and dwarf 
spheroidal galaxies at galactocentric distances from 20 up to 
60 kpc has been analysed by Hartwick and Sargent (1978). 
They determine the mass within 60 kpc, under the two oppo- 
site extreme assumptions of purely radial or purely isotropic 
velocity distributions, respectively. Translated into a circular 
velocity, this leads respectively to w rot = 159 ± 30km. s _1 
and v rot = 237 ± 30km.s _1 . Caldwell and Ostriker (1981) 
used the more realistic geometrical mean of the two values 
^rot = 192 ± 42 km.s -1 . The latter finally has been scaled by 
Bahcall et al (1983) to w rot (60kpc) = 205 km.s -1 (to which 
must be added a ±45 km.s -1 uncertainty) in order to take into 
account more recent determinations of the local parameters 
near the Sun. These results are consistent with the more re- 
cent and much more precise (although requiring a modelling of 
the galaxy profiles) determination of the mass of the Galaxy by 
Kochanek (1996), i.e. 3.8 x 10 11 M < M 30 < 5.4 x 10 11 M 
and 4.0 x 10 11 M < M 60 < 8.7 x 10 11 M at the 90% C.L., 
where M30 and Mqq denote the mass within 30 and 60 kpc, re- 
spectively. This yields, within the la error level we have been 
using all along the present analysis, v mt (30) = 257±8 km.s -1 
and u ro t(60) = 210 ± 13 km.s -1 . 

Fich and Tremaine (1991) discuss the implications of the 
velocities measured for the "Magellanic Stream", a hydrogen 
bridge between the Magellanic Clouds and the Galaxy. A de- 
tailed modelling shows that the r -2 (r s fits better) halo 
should extend up to at least 60 kpc. Although such a conclu- 
sion is reached by means of several "natural" assumptions, and 
may be questioned, the raw data show gas velocities around 
200 km.s -1 . 

Peebles (1989, 1990, 1994) has examined the total mass 
of the Local Group, dominated by Andromeda and the Milky 
Way. He finds total masses around 4 x 10 12 M , suggesting a 
total mass of about 2 x 10 12 M for our Galaxy. Moreover, a 
modelling of all objects at a distance of 4 Mpc from the Galaxy 
shows indeed that their velocity distribution is consistent with 
the mass being attached to the objects. 



In summary, all these studies lead to the conclusion that the 
Milky Way has a mass ~ 2 x 10 12 M , with a nearly con- 
stant rotation velocity t> rot = 220 ± 20 km.s -1 , which implies 
an extension of a Galactic dark halo up to at least 100 kpc. This 
undetected mass has to be compared with the mass under the 
form of visible matter, ~ 10 11 M , which sets the scale for the 
amount of dark matter in the Galaxy. 



velocity dispersion (v z ) and density of one tracer population 
(Binney & Tremaine, 1987) by: 



K z {z) = -- dp{vl) 



p dz 



(1) 



By Taylor-expanding the Galactic potential above the disk, 
<&(z) w K z z + Fz 2 , the Poisson-Jeans equation yields (Binney 
& Tremaine 1987): 



<9$ 

K z (z) = =K zQ + 2zF 
oz 



(2) 



where the coefficient K z0 depends on the mass content in the 
disk (assumed to be infinitely thin, which means the above ex- 
pression is good for z ^ 300 pc only) while F — 2-KGphaio is 
the halo contribution. The local surface density S in the disk 
is related to K z q by: 



So = 



K zQ 
2ttG 



(3) 



The surface density near the Galactic mid-plane determines 
the mass of the disk, for a given disk scale length R^. 



M disk «2^£ i? 2 e fl °/^ 



(4) 



In practice, K z (z) is measured at some height z and must 
be corrected from the halo contribution 2zF to determine K z q 
and E disk- Bahcall (1984a, 1984b), in agreement with Oort 
(1960) found S ps 7OM pc -2 with a claimed small er- 
ror. Bienayme, Robin and Creze (1987), using a sampling ex- 
tending to higher latitudes, obtained a better determination of 
K z (z) from 100 pc to 1 kpc, based on a modelling of the disk 
and the halo stellar population arising from available observa- 
tions. They found a lower surface density for the disk, 



50 M pc" 



(5) 



with a 10 M pc -2 uncertainty, and stressed that the uncertain- 
ties in Bahcall's extrapolation to low z were underestimated. 
These authors, however, actually measure the acceleration at 1 
kpc, where the surface density is essentially the total surface 
density, and get K z (l kpc)/27rG « 7OM pc -2 . The afore- 
value of S quoted for the disk is obtained after substracting 
the standard halo contribution^ (see eqn.(2)), and thus becomes 
model-dependent. Subsequently, Kuijken and Gilmore (1989, 
1991), found an observed total density £(1.1 kpc) = 71 ± 6 
M pc -2 and argued that, after a theoretical correction for 
the contribution of a standard massive halo similar to the one 
used by Bienayme et al. (1987), the surface density of the disk 
should be reduced to : 



2.2. Local dynamical surface density 

An important constraint on the mass density of the Galactic 
disk in the solar neighbourhood is obtained from the study of 
the vertical acceleration K z by means of the local stellar ve- 
locity distribution. K z (z) is related to the (measured) vertical 



S = 48 ± 9M pc" 



(6) 



1 Throughout this paper, the standard halo is defined as a halo with 
a density profile p(r) — r 2 +R i > wrtn v TOt — 220 km.s -1 and 
R c — 5 kpc, see paper II. Note that this halo is slightly heavier than 
the one used for instance by the MACHO group, which corresponds 
to v TOt = 204 km.s -1 . 
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in remarkable agreement with the afore-mentioned determina- 
tion. A subsequent reanalysis of the Kuijken & Gilmore (1989) 
work gives in fact S = 54 ± 8M pc" 2 (Gould, 1990). We 
stress, however, that the two afore-mentioned surface densities 
at 300 pc depend somehow on the assumption made for the 
dark halo contribution (Eq. (fy). 

More recently, Bahcall, Flynn and Gould (1992), using the 
new data of Flynn and Freeman (1993), re-analyzed Bahcall's 
(1984) previous determination, and obtained 

S Q = 85±25M pc~ 2 . (7) 

The new lower limit weakens substantially the conflict with 
the afore-mentioned values, at the price however of a large un- 
certainty. The most recent analysis of these data by Flynn and 
Fuchs (1994), who added a new normalization point, yields a 
best fit model: 

S Q = 52 ± 13 M Q pc- 2 (8) 

Both groups include a standard halo contribution in their deter- 
mination but for the average height of their sample (~ 300 pc), 
the halo correction is small (~ 6M pc~ 2 ) and the associated 
uncertainty lies well within the observational error bars. The 
reason why Bahcall et al. (1992) considered their value to dis- 
agree significantly with the Bienayme et al. (1987) and Kuijken 
& Gilmore (1991) values is that even with a la deviation, the 
probability for S , as given by (7), to be less than 60 M pc -2 
is only 7% and thus has only one chance in 14 to occur, as 
stated by these authors. More optimistically, comforted by the 
recent Flynn & Fuchs (1994) result, we note that all these de- 
terminations of S should rather be considered as consistent 
at the lcr level: the uncertainty-weighted average of these four 
different results (5)-(8) yields: 

S = 51±6M pc- 2 , (9) 

very close to Eq. We thus may consider the value as 
the most accurate present determination of the solar dynamical 
density. 

3. Star counts 

3.1. Disk 

The determination of the faint end of the stellar luminosity 
function (LF), an essential issue to infer the low-mass star and 
the sub-stellar contributions to the Galactic mass budget, has 
also been subject to strong debate. Wielen et al. (1983) found 
that, whereas the LF was steadily rising up to My = 13, there 
was a dip beyond this value, leaving very little possibility for 
a substantial contribution of low-mass stars to the mass of the 
disk. Subsequent studies (Stobie et al. 1989, Tinney et al. 1992) 

2 Strickly speaking, the Bahcall et al. (1992) and the Flynn & Fuchs 
(1994) values are not completely independent but the weight of the 
Bahcall et al. value is small. In any event, the value used as a reference 
all along our calculations is the most recent Flynn & Fuchs one. 



showed that the dip is even more pronounced and starts already 
at My = 12. Gould, Bahcall and Flynn (1996;1997), using 
HST observations at faint magnitudes, confirm a dropping LF 
beyond My = 12, which corresponds to m ~ 0.25 M , al- 
though the last bin of the HST LF clearly shows an upturn at 
My > 16. Most of these LFs are based on a photometric deter- 
mination of the distance, which thus relies on a color-absolute 
magnitude relationship, up to about a few hundred pc. On the 
other hand, Kroupa (1995) used a nearby LF, allowing a geo- 
metric determination of the distance. This yields larger statis- 
tical errors, due to the number-limited sample, but better sys- 
tematic corrections, in particular the ones arising from the un- 
resolved binary companions. The nearby LF agrees with the 
photometric one up to My ^12 but is essentially flat beyond 
this limit. Doing a careful correction for the Malmquist bias 
and for the presence of unresolved binaries, Kroupa, Tout and 
Gilmore (1990; 1993) and Kroupa (1995) have shown that the 
properly corrected ground-based photometric LFs are consis- 
tent with the local geometric LF. 

The issue, however, remains unsettled for the HST LF 
(Gould et al., 1996; 1997). The HST is almost not subject to 
Malmquist bias since its limit magnitude allows the observa- 
tion of the bottom of the main sequence up to the edge of the 
disk, but it is subject to uncompleteness due to unresolved sys- 
tems near the faint end[| Gould et al. (1997), however, found 
that even with this correction, their LF is inconsistent with the 
nearby LF. Near the end of the LF, the binary correction is, ac- 
cording to these authors, at most a factor 2, whereas a factor 
~ 5 is required to reconcile the two kinds of LFs. Therefore, 
unresolved binaries are certainly not the source of the discrep- 
ancy between the two LF's. An important quantity to calculate 
correctly the correction due to binaries is the rate of binaries 
and the mass-distribution of the secondaries. An accurate deter- 
mination of such quantities require decade-long observations, 
as conducted for example by Duquennoy and Mayor (1991; 
DM91) for nearby F and G stars, with a fairly accurate sensitiv- 
ity (~ 300 m.s -1 ). Note that, altough there are severe natural 
observational biases toward the detection of near equal-mass 
systems either from imaging or from spectroscopic surveys, the 
completeness and the uniformity of the DM91 sample render it 
almost without any observational selection effect, yielding an 
almost unbiased mass-ratio determination for F and G star bi- 
naries. The main conclusion of the DM91 study is that a fairly 
large number (~ 60 — 70%) of single stars form a multiple 
system, and that the shape of the mass-distribution increases 
toward small mass ratios (q = m,2/mi < 1; see their figure 
10). 

M-dwarf surveys to address the same issue are still in their 
infancy and suffer from severe incompleteness and observa- 
tional biases, so that it will require more years to do a simi- 
lar analysis for M-dwarf binaries. In any case, as mentioned 
above, even a full knowledge of the binary correction will not 
account for all the discrepancy between photometric and geo- 

3 The HST is insensitive to binaries with separations ^ 0. 3 and 
thus misses essentially all secondaries in late M-dwarf binary systems. 
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Fig. 1. Mass functions derived from Kroupa (1995; solid line) 
and Gould et al. (1997; dashed line) LFs. Also shown are the 
fits for both MFs, which are identical for m > 0.35M Q (dotted 
lines). See table |l]for the parameters of the fits. The dashed line 
is an overall fit given by Eq. ( |To| ) with a = 2, which reproduces 
reasonably well the MF derived from Kroupa's LF (note that 
there are only 2 stars in the lowest bin of Kroupa's LF; see 
text). 

metric LFs. For this reason, we have considered in the present 
paper both (nearby and HST) LFs. 

We have converted both LFs into mass-functions (MF), us- 
ing low-mass star models which accurately reproduce observed 
mass-My relationships (Chabrier, Baraffe & Plez, 1996). Both 
MFs are displayed in Figure [I] in a log-log scale. The afore- 
mentioned discrepancy between the two LFs is obvious on the 
MF, in particular in the domain ~ 0.1 — 0.25 M Q , with a factor 
~ 5 ratio at ~ 0.1 M . 

Noting that the two lowest bins in the nearby LFs are the 
least statistically significant (2 stars, see Kroupa (1995)), Mera 
et al. (1996a) parametrized the disk low-mass star MF in the so- 
lar neighborhood down to the vicinity of the hydrogen-burning 
limit by the following form: 

M (m) « 1.5 ± 0.4(^-)- 2±0 - 5 Mq 1 pc" 3 , (10) 

The undetermination of the exponent reflects the fact that 
this fit is a reasonnable overall power-law representation of the 
true MF but is not perfect. As shown by Kroupa et al. (1993), 
the MF determined from the nearby LF is better described by 
a series of segment power-law fonctions /i(m) = dN /dm cx 
mT a . We then have fitted the MFs derived from the nearby and 
the HST LFs with such segmented power-laws. The character- 
istics of these MFs are given in table [j] and illustrated in Figure 
1. We will refer to these fits as MF (a) for Kroupa (1995) and 
MF (b) for Gould et al. (1997). 

If we exclude the last bin, the HST MF corresponds to a 
slope a < 1.0 for m/M Q 0.6. However, an extrapolation 
of the MF (ignoring the last bin) all the way below 0.2 M© 
would correspond to about 2 stars in the last bin whereas 10 



Table 1. Fit of the two MFs considered in this paper. The MF 

(a) has been derived from Kroupa (1995) LF, whereas the MF 

(b) has been derived from Gould et al. (1997) LF. The MFs and 
their fits are displayed in figure [l] Each MF is considered to 
be a 3-segment power-law /i(m) = dN/dm = AmT a . Both 
MFs are identical for m > 0.35 M©. The mass is in M©, the 
mass function in Mq 1 pc -3 . The upper mass limit has been set 
arbitrarily to 10 M©, since it has almost no influence on the 
mass density. Note that 0.4 should be added to the slope de- 
rived from the HST LF for m < 0.6 M© to account for missed 
binaries (Gould et al., 1997) 
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are seen. The probability to see 10 stars whereas 2 are expected 
is less than 5 x 10~ 5 . This suggests that the last bin in the 
HST LF is indeed significant, adding credibility for a rising MF 
near the BD limit. As shown in Figure [j], the fit does not repro- 
duce exactly the last bin. This takes into account a possible bias 
which tends to overestimate the last bin and to underestimate 
the penultimate one (Gould et al. 1997). 

Indeed, for both MF's (nearby and HST), there is a hint for 
a rise near 0. 1 Msol, suggesting the possibility for a substan- 
tial amount of brown dwarfs in the disk. We will take this as a 
basis for extrapolating the stellar MF into the brown dwarf do- 
main. This is consistent with a recent study by Mazeh, Latham 
& Stefanik (1996) on the distribution of secondary masses near 
the substellar limit. Although this study relies on the detection 
of 3 very-low-mass companions out of 20 F and G dwarfs and 
subgiants, and thus preclude the precise determination of a sub- 
stellar MF, these detections suggest that the mass distribution 
(the primaries of the sample all have ~ 1M Q ) rises near the 
sub-stellar limit. 

In any case, whereas the precise determination of the shape 
of the MF near the bottom of the main sequence, which requires 
larger statistics, is of prime importance for star-formation the- 
ory, the qualitative behaviour, and its overall characterization, 
is what is most essential to estimate the stellar and sub-stellar 
mass contributions to the disk mass budget. Both rising MF's 
near the BD limit suggest a substantial amount of substellar 
objects and thus a contribution to the disk mass density. The 
extension into the brown dwarf domain will be analyzed in the 
next section, in connection with microlensing observations. 

The integration of the nearby MF displayed in Fig. [j] gives 
the local low mass star (0.07M© < m < 0.65M©) densi- 
ties plms = 2.56 ± 0.25 x 10 -2 M©pc -3 in case (a), and 
p LMS = 1.43 ± 0.09 x 10" 2 M©pc" 3 in case (b). Since 
the HST LF does not include companions of multiple systems, 
we have tentatively corrected for this bias on the basis of the 
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Duquennoy and Mayor (1991) observed mass ratio distribution 
mentionned above. The correction is approximated by a a lin- 
ear fonction, whose value is 10% at 0.6 M , and 80% at 0.1 
Mq. We renormalize the HST MF to keep the same normal- 
ization at 0.6 M Q with this correction. The resulting density is 
Plms — 1-88 x 10~ 2 M Q pc~ 3 , almost in agreement with the 
nearby value. An average value is then: 

Plms = 2.2 ± 0.3 x 10" 2 M pc~ 3 

With a vertical sech 2 density profile of scale height hd ~ 
320 pc (Gould et al., 1997), the contribution of low-mass stars 
(LMS) down to the bottom of the main sequence (0.1 ^ 
m/M < 0.6) to the disk surface density is then: 

E lms «14 ± 2M pc- 2 (11) 

Adding the contribution of more massive stars, obtained with 
the MF and the luminosity-dependent scale height determined 
by Miller & Scalo (1979), converted into a mass-dependent 
scale height, yields the total contribution of main sequence 
stars: 

E* = 24 ± 3M Q pc~ 2 , (12) 

Including the stellar remnants, white dwarfs (Liebert, Dahn 
and Monet, 1988) and neutron stars, E w d+ n s = 2 to 4 
M pc~ 2 , and an estimated interstellar gas contribution E gas = 
10 to 15 M Q pc~ 2 (Bahcall, 1984), the total surface mass den- 
sity under the form of directly detected baryonic components 
in the disk is thus: 

(S ) vis =40 ± 4M Q pc- 2 (13) 

in good agreement with Gould et al. (1997). The limits are de- 
rived from very conservative estimates for the errors. Using a 
sech(|z|/zo) 2 vertical density distribution instead of an expo- 
nential profile leaves these values almost unchanged. 

3.2. Bulge 

The mass and the shape of the Galactic bulge will be discussed 
in more details in the next section. This mass is estimated to 
be ~ 2 x 10 10 M Q , from dynamical considerations (Zhao et 
al., 1996), with a size of about 1 kpc (Kent 1992). The faint 
end of the LF in the Baade's Window, down to My ~ 10, 
i.e. m « 0.5 Mq, has been studied with the HST Wide Field 
Camera (WFC) by Holtzman et al (1993). These authors find 
a Salpeter slope for the MF x w 2.35 ± 1, over the afore- 
mentioned mass range. There is no indication for this MF to 
be different, except for its local normalization, from the one 
determined for the disk, Eq. (|l0|). 

3.3. Thick disk 

A disk+bulge model does not reproduce successfully the ob- 
served star counts and kinematic properties. This led Gilmore 
& Reid (1983) to introduce the so-called thick disk, with a scale 



height /ltd ~ 1 — 1.5 kpc, which seems to be necessary for an 
accurate galactic modelling. Since the stellar velocity disper- 
sion {v 2 z ) and the oblateness of the disk are related, a thick disk, 
by definition, is less elongated than the young (thin) disk, and 
has a larger velocity dispersion. This yields a better agreement 
with observed velocity dispersions, and with star counts in gen- 
eral. As shown recently by Reid et al. (1996) and Chabrier & 
Mera (1997), a thick disk with /i t d ~ 1 kpc and q w 0.6 is con- 
sistent with the Hubble Deep Field counts at faint magnitudes. 
A larger scale height and a spherically symmetric distribution 
predict substantially larger star counts. In spite of the small vol- 
ume density of the thick disk in the solar neighborhood, about 
a few percents of the young disk at 1 kpc (Gilmore, Wyse & 
Kuijken; 1989), its contribution to the dynamics is not negli- 

/+oo 
Ptd(R, z)dz 
-oo 

represents ~ 10-20% of the young disk contribution. 

Taking into account the thick disk contribution, and using 
the parameters derived by Gould et al. (1997), the disk stellar 
surface density is raised by ~ 3 M pc~ 2 to 

E* = 27 ± 4M pc- 2 , (14) 

so that the detected surface density ( pj| ) reads: 

(E ) vis = 43 ±5M pc- 2 (15) 

3.4. Spheroid 

The spheroid is also called in the literature the stellar halo. 
It is defined by a spherical density distribution decreasing as 
p{r) oc r~ 3 (Hubble profile) or equivalently a r 1 / 4 , de Vau- 
couleurs profile. It differs from the central bulge essentially 
by a lower metallicity, [M/H] w -1.5 (Monet et al., 1992; 
Leggett, 1992; Baraffe et al., 1995, 1997) and larger velocity 
dispersions, ctr w 160 km.s -1 , a z = ag w 90 km.s -1 (Dahn 
et al. 1995; Beers & Sommer-Larsen, 1995). The mass of the 
spheroid is ~ 10 9 M (Bahcall, 1986) so its contribution to the 
galactic mass is not important. Its main contribution concerns 
the halo star counts. 

Richer and Fahlman (1992) obtained for the spheroid a 
steeply rising LF down to the end of the observations whereas 
Dahn et al. (1995), from observations of high-velocity («tan > 
200 km.s -1 ) subdwarfs in the solar neighorhood, get a LF 
which decreases at the faint end (My ^12). 

Mera, Chabrier & Schaeffer (1996) and Chabrier & Mera 
(1997) (see also Graff & Freese, 1996) have determined the 
MF of the spheroid from these observed LFs and from mass- 
magnitude relationships derived from LMS models which re- 
produce accurately the main sequences of metal-poor globu- 
lar clusters, observed with the HST deep photometry surveys, 
down to the vicinity of the hydrogen burning limit (Baraffe et 
al., 1997). The spheroid MF in the solar neighborhood deduced 
from the Dahn et al. (1995) subdwarf LF is reasonably well 
parametrized as (Chabrier & Mera, 1997): 

/ \ -1.7±0.2 

M (m) = 4.0 ± 1.0 10- 3 iyoY^J M^ 1 pc" 3 , (16) 
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Integration of this MF gives the halo main sequence (to < 
0.8 M Q for t > 10 Gyr) stellar density in the solar neighbor- 
hood : 

p K re 1.0 x 10" 4 M Q pc- 3 (17) 

As shown by Chabrier & Mera (1997), star count predic- 
tions based on this MF are in excellent agreement with the 
Hubble Deep Field observations at very faint magnitudes both 
in the V and /-bands. Comparison of this value with the disk 
main-sequence stellar density determined in §3. 1 gives the nor- 
malization of the spheroid/disk stellar population at 0.1 M Q , 
namely p* sph /p* d ~ 1/400 (Chabrier & Mera, 1997). 

3.5. Dark Halo 

Throughout this paper, the term "dark halo" defines what is 
characterized by the p(r) oc r~ 2 density-profile, a consequence 
of the flat rotation curve in the outer part of the Galaxy (see 
§2). The existence of a halo in the Galactic structure is ren- 
dered necessary i) to account for the total mass of the Galaxy 
(see §2.1.2), since the disk(s), bulge and spheroid total mass 
amounts at most to ~ 10 11 M© and ii) for stability conditions 
(Ostriker and Peebles, 1973), although the presence of a cen- 
tral bar weakens this latter argument. There is observational 
evidence that galactic dark halos present some oblateness with 
q ~ 0.6 (Sackett et al. 1994). Were the mass-distribution of the 
dark halo 2-dimensional (i.e. under the form of a disk), the sur- 
face density in the solar neighborhood, E Q pa = 
210 Mq pc~ 2 would be a factor 3 larger than the dynamical up- 
per limit (see §2.2) (see Mera et al., 1997, paper II). This con- 
dition, and the motion of galaxies and gas in the Local Group 
(§2.1) are strong arguments in favor of a large, 3-dimensional 
halo around the centre of the Galaxy. 

The MF in the halo is presently unknown. The observation 
of the bottom of the stellar main sequence, m ~ 0.1 M Q , at a 
distance of 20 kpc, about the limit of the spheroid, requires sur- 
veys down to apparent magnitudes / > 28 (Baraffe et al. 1997). 
It requires also an excellent angular resolution (< 0.1 arcsec) 
to distinguish stars from galaxies. This is the reason why early 
attempts using ground-based observations were never brought 
to a conclusive end. The recent Hubble Deep Field star counts 
(Flynn et al., 1996; Reid et al., 1996; Mendez et al., 1996) up 
to / = 28 are insufficient to determine a halo LF, because of 
the limited field of view, but represent a severe constraint on 
the Galactic model. These faint magnitude star counts have 
been shown to be entirely consistent with a thick disk + flat- 
tened spheroid population, characterized by their respective I- 
(V-I) relationships and mass functions ( JTo| ) and ( p"rj| ) (Chabrier 
& Mera, 1997). As shown by these authors, a (even flattened) 
dark halo with the afore-mentioned 1/r 2 profile and MF (16) 
would predict at least 500 more stars than observed in the HDF 
field of view. This shows that, if the dark halo LF is the same 
as the spheroid one, there is a negligible stellar population in 
the dark halo, less than 0.1% of the dynamical mass. 

A promising technique to determine more accurately the 
halo MF comes from the recent surveys of the outer parts of 



nearby galaxies (Sackett et al. 1994; Lequeux et al. 1996) and 
from the relation between the slope of the MF and the inferred 
colors of the halo diffuse emission (Mera, 1997). 

Another important indication is given by the LFs of glob- 
ular clusters. The MFs of the different clusters observed with 
the HST have been determined recently by Chabrier & Mera 
(1997). The typical LFs, determined near the respective half- 
mass radii, are rising up to Mboi ~ 10 and drop below. This is 
consistent with monotonically rising MFs down to ~ 0.1 Mq, 
with a very weak dependence on metallicity, well described by 
power-law functions with a ~ 0.5 — 1.5, similar to the ones de- 
termined in the previous sections for the disk and the spheroid. 
Although mass segregation effects are expected to be small 
near the half-mass radius (see e.g. King et al., 1995; Chabrier 
& Mera, 1997), the evolution of the clusters, and their evapo- 
ration, might lead to initial MFs (IMF) slightly steeper that the 
afore-mentioned determination. This would yield even better 
agreement with the MF (|l6|) and would be consistent with the 
suggestion that the formation of spheroid field stars and glob- 
ular clusters has occured in a similar manner and time frame 
(Fall and Rees, 1985). 

Since the dark-halo stellar population is likely to be in- 
significant, as mentioned above, Eq. (|l6|) gives the total halo 
(spheroid + dark halo) stellar contribution to the dynamically- 
determined local density pdyn ~ O.OlM .pc -3 , namely 
Ph+ /pdyn ~ 1%. This puts more stringent limits than the ones 
obtained from the HST counts alone (Flynn et al., 1996) and 
shows convincingly that low-mass stars represent a completely 
negligible fraction of the halo mass. 

In summary, the star count analysis of different stellar pop- 
ulations corresponding to different regions of the Galaxy, based 
on accurate evolutionary models for low-mass stars (Baraffe et 
al., 1997; Chabrier et al., 1996, 1997), leads to the determina- 
tion, slope and normalization, of the MF in the disk, the bulge 
and the spheroid, down to the bottom of the main sequence. 
Surprinsingly enough, given the different metallicities in these 
regions and the wide range of stellar masses considered, these 
MF's are very similar and reasonably well described by power- 
law MFs dN/dm oc m~ a with a ps 1.5 - 2. A flatter MF be- 
low ^0.5 M for the disk, as sugested by the HST LF, cannot 
be excluded but a more precise determination requires better 
statistics at the faint end of the LF. These MFs, consistent with 
all observed LFs and deep-photometry counts, yield a fairly 
well determined contribution of all main sequence stars to the 
disk and halo mass budgets, the dark halo stellar contribution 
to the Galactic mass being essentially insignificant. The deter- 
mination of the amount of mass under the form of sub-stellar 
objects requires the analysis of the microlensing experiments. 

4. Microlensing 

Within the past few years the microlensing technique has been 
applied succesfully to the search for dark matter, by infer- 
ring the presence of dark objects in the Galaxy through their 
gravitational effect on luminous matter, as proposed initially 
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by Paczynski (1986). Several groups worldwide have carried 
out microlensing experiments to determine the amount of dark 
matter either in the halo of the Galaxy (EROS and MACHO 
experiments) or in its central parts (OGLE and MACHO ex- 
periments). We examine below these various experiments and 
connect the results with the afore-mentioned stellar MFs to de- 
termine the contribution of sub-stellar objects to the galactic 
mass budget. 

4.1. Disk and Bulge 

Numerous events have been detected in the direction of the 
Galactic centre by the OGLE (Paczynski et al., 1996) and MA- 
CHO (Alcock et al. 1997) surveys. We focus here on the 40 
events^ obtained from the first-year analysis of the MACHO 
collaboration, whose detection efficiency e(t) is well deter- 
mined (Alcock et al., 1997). The number of events observed 
by OGLE during the first two years is too small and the at- 
tempt to use averages would be dominated by small number 
statistics. At any rate, adding, as an extreme value, the 9 OGLE 
events to the 40 MACHO ones would not improve significantly 
the statistics. 

The interpretation of the three large-time events (t e > 75 
days) in the MACHO survey is still unclear and will be dis- 
cussed below. Note that there is no event with 35 < t e < 75 
days. 

The number of observed events allows a statistical analysis 
from their time distribution. The average effective time (t e ) can 
be used to estimate the average lens mass (to) and thus the 
minimum mass down to which the mass function ( fic| ) extends 
(see Appendix): 



(V^) = — =<o<— r^v/^o^r 1 (is) 

where the averages are relative to the effective probability dis- 
tributions Pcff for each variable (see Appendix). For a simple 
gaussian velocity distribution, = cry/w/2, where a is 

the velocity dispersion. The geometric factor (\J x(\ — x)) is 
approximatly 0.39 for a standard (Bahcall & Soneira, 1980 + 
Kent, 1992) disk+bulge model, and is not very sensitive to the 
choice of the model. The mean duration of the 40 events, tak- 
ing into account the experimental efficiency, is (t e ) = 16.8 
days (omiting the efficiency leads to (t e ) = 20 days). Then a 
first estimate of the average lens mass is: 

(V^) 2 = 0.15M Q (19) 

which is only a rough estimate, with a velocity dispersion 
a = lOOkm.s -1 for the bulge. Excluding the three longest 
events yields (\/m) 2 = 0.09M©. These estimates introduce 
the dilemma that will be discussed in detail below: in the first 
case, microlensing results are compatible with essentially no 

4 We have excluded the two events that fail MACHO cuts (104-B 
and 1 1 1-B), two events suspected to be variable stars (1 13-C and 121- 
B) and the binary event which is not accounted for by the efficiency. 



brown dwarf in the disk/bulge, whereas in the second case, the 
mass function must extend significantly in the brown dwarf do- 
main. 

As mentioned in the Appendix A. 2, we stress that P(t e ) 
has a long power-law tail at large times, since t e oc wj 1 , which 
yields an unphysical divergence of the moments. Therefore, 
the statistics for the mass of the microlensing events can not 
be evaluated correctly from the calculation of the moments, as 
done sometimes in the literature, but must be evaluated accu- 
rately by using the detailed time distribution (Note, however, 
that the moments of 1 /t e behave correctly). 

We have calculated the time distribution with a Monte- 
Carlo algorithm (see Appendix), which determines the number 
of events: 

N th = Ex / e (t e )—dt e (20) 

Jo Me 

where E is the total exposure (in star-years) and e the exper- 
imental efficiency factor (Alcock et al., 1997). Calculations 
have been done for a standard disk model (Bahcall & Soneira, 
1980), and a Kent (1992) model for the bulge, for the two 
MFs considered in §3.1. We have taken velocity dispersions 
a ~ 100 km.s -1 for the bulge and a ~ 20 km.s -1 for the disk 
(40 km.s -1 in the radial direction), respectively, with a linear 
interpolation in the intermediate region located between ~ 2 
and 3 kpc from the Galactic center. As will be shown below, 
the uncertainty in the velocity distribution of the stars is the 
main cause of uncertainty in the present calculations. 

Our calculations take into account the motion of the Sun 
and of the source star in the determination of the lens veloc- 
ity as well as the variation of the distance of the source stars 
in the disk and the bulge (see Appendix). The distance of the 
source stars in the bulge is typically 7.5 kpc from the Sun, as 
determined by Holtzman et al. (1993) from the observed LF in 
Baade's Window. The time-distribution and number of events 
obtained from these calculations, for a given model, are com- 
pared to the MACHO observations (Alcock et al., 1997). Since 
the efficiency does not include binary events, we must exclude 
those events. However, only a small fraction of the events are 
expected to be due to both components of a binary system (1 
out of 45 in the first year, a fraction which seems to be con- 
firmed by the following two years). In most cases, the Einstein 
disks of the stars in a multiple system do not overlap, and can 
in first approximation be treated independently. The time distri- 
bution is therefore barely affected by the binarity. In any case 
the number of events has to be corrected only slightly (typi- 
cally a few percent). In this section, we are mostly interested in 
the determination of the mass function in the substellar regime, 
which is derived only from the time distribution of the observed 
events. The consistency between models and observations is 
analyzed with a Kolmogorov-Smirnov (KS) testQ. 

The KS test does not give a confidence level, but rather an ex- 
clusion criterium. Its significance function has a very rapid variation 
between 10% (distributions most probably different) and 90% (distri- 
butions most probably identical), and any value in between does not 
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We will first assume that the mass function in the bulge 
is the same as in the solar neighborhood. The mass function 
is extrapolated in the brown dwarf domain down to a mini- 
mum mass, mj n f, which is determined when the best agree- 
ment between theory and observation is achieved. The results 
are shown in Figure ^, for the mass functions (a) and (b) defined 
in §3.1, with minf = 0.075 M Q and m m f — 0.05 M Q , respec- 
tively. The KS test gives a probability of 85% for MF (a) and 
only 45% for MF (b). The value m inf = 0.05 M for MF (b), 
and its steeper slope in the brown dwarf domain renders MF (b) 
very similar to a Dirac delta-function. Longer timescale events 
are then more difficult to reproduce. On the other hand, for MF 

(a) , the relative contribution of higher masses is larger, with a 
better fit of the long timescale observations. 

We have also computed the time distribution in case (b) 
with mjnf = 0. IMq, which correspond to the absence of brown 
dwarfs in the disk. This corresponds to the general decreasing 
behaviour of the HST MF, ignoring the last bin. The time distri- 
bution, also shown in Fig. |2| is incompatible with the observed 
time distribution, with a KS result of 5%, as already pointed 
out by Han and Gould (1995). These results suggest that brown 
dwarfs are required to explain all the 40 events of the MACHO 
first year bulge observations, unless the effect of amplification 
bias (Alard, 1997; Han, 1997) is important enough to account 
for the excess of short timescale events. 

It has been suggested (Alcock et al 1997; Han and Gould, 
1995) that the large time events may have another origin than 
main sequence stars, and might be due to stellar remnants, 
white dwarfs or neutron stars. They may be due also to rare 
statistical fluctuations due to the large tail of the time distri- 
bution or, more speculatively, to dwarf nova eruptions (Delia 
Valle and Livio, 1996). If we exclude the 3 longest events from 
the present analysis, the agreement between theory and obser- 
vation improves significantly, as shown in Fig. || for the same 
MFs. This yields minimal masses for the two mass functions 
TOi n f = 0.056 Mq in case (a), and minf = 0.047 M in case 

(b) . The KS significance of the two MFs is 50% in case (a) and 
40% in case (b). 

The minimum masses given above depend on the model 
used for the calculations, in particular the poorly-determined 
velocity dispersions in the bulge population (Ibata & Gilmore, 
1995). This yields the largest uncertainties in the present re- 
sults, since (m) oc (a 2 ). 

These mass functions yield a brown dwarf contribution to 
the surface density in the solar neighborhood, assuming the 
same scale height as for the M-dwarf population, i.e. hd ~ 300 
pc: 



Sbd ~ 3 ± 3 M pc" 



(21) 



Therefore the brown dwarf contribution to the local surface 
density is about 10% to 20% of the total stellar contribution 
(14). This represents ~ 10% of the observed baryonic mass 




Fig. 2. Observed time distribution (histograms) for all MACHO 
events compared with the distribution of two models. Solid 
line: MF (a) with m inf = 0.075 M Q , dashed line: MF (b) with 
m inf = 0.05 M , and dotted line: MF(b) with m inf = 0.1 M . 
The latter case is excluded by a KS test (see text), showing that 
brown dwarfs are required to explain the observed distribution. 




Fig. 3. Same as figure g, without the three longest events (see 
text). Models are: MF (a) with m in{ = 0.056 M Q (solid line), 
MF (b) with m in £ = 0.047 M Q (dashed line). 



(|l5]). As shown in Eq. (pi]), the brown dwarf contribution is 
also compatible with zero. The total surface density in the form 
of baryonic matter in the disk is thus (Eqs. ( |T5| ) and (pT|)): 



-'baryon 



46 ± 6M pc~ 



(22) 



give a strong conclusion. But the KS test is better than the \ test for 
low statistics, and does not depend on any binning of the data. 



consistent with the dynamical determination (|8J). There is no 
need for any additional dark matter. The contribution of stellar 
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and substellar objects to the surface density is (Eq.(|l4]) and 
Eq.@): 



(23) 



E ss = 30 ± 5 M pc 



This yields an optical depth (see e.g. Kiraga & Paczyhski, 
1994): 



2ttG 



Tdisk 



3c 



2 p*i?Q 



5.5xlO^ Ess + 5OM0 f 
50M Q pc- 2 

9.0 ±0.5 x 10" 7 , 



(24) 



i.e. Tdi S k ^ 10 6 , a factor ~ 2 to 3 smaller than the observed 
value T obs ~ 2.4 ± 0.5 x 10~ 6 (Alcocket al., 1997). Therefore, 
although the observed time-distribution of the events is well 
reproduced by the present model, the optical depth amounts 
to ~ 1/3 of the observed value. This clearly suggests a bulge 
strongly elongated along the line-of-sight or, similarly, a bar 
model (Zhao 1994). We have not included a bar-model in our 
calculations, but from the present analysis, its contribution to 
the optical depth corresponds to: 



T ba 



1 x 10- 6 M bar /(10 10 M Q ) 



(25) 



in good agreement with various models (Zhao, Spergel & Rich, 
1995, 1996; Han & Gould, 1995; Stanek et al., 1997; Bis- 
santz et al., 1997). Infrared observations from DIRBE seem 
to constrain the mass of the bulge within 2.4 kpc to M& ~ 
7.2— 8.6x 10 9 Mq, although this value is obtained by substract- 
ing a model-disk contribution; a more general constraint con- 
cerns the bulge+disk mass M b+d (< 2.4 kpc) = 1.9 x 10 10 M o 
(Bissantz et al., 1996). Note however that the three long-time 
events contribute about 1/3 to the observed optical depth, so 
that the discrepancy between the value ( j24| ) and the observed 
one is reduced appreciably if these events are removed, leading 
to a smaller mass for the bulge. Moreover, as noted recently 
by Han (1997) and Alard (1997), the optical depth toward the 
bulge might be overestimated by as much as a factor 1 .7, be- 
cause of amplification bias due to pixel lensing. 



4.1.1. Events whose source is a bulge clump giant 

As mentioned by Alcock et al. (1997), there are significant un- 
certainties in the estimated bulge optical depth, because essen- 
tially of blending effects. For this reason, the MACHO collab- 
oration has conducted a separate analysis of the events whose 
source is a bulge giant (Alcock et al. 1997). Although the statis- 
tics are significantly diminished (1.3 x 10 6 source giants and 
13 events), the efficiency and the reliability of the results are 
improved appreciably. The corresponding optical depth is also 
substantially larger, but only at the ~ ler level, than the value 
obtained when including the totality of the events, namely 



' giant 



(te 



3.9 x 10 6 . The mean duration of these events is 
33.8 days, which yields an unrealistic average mass of 



0.6 Mq from relation ([TJ]). This high value is confirmed by the 
complete calculation, with a minimum mass mj n f » 0.3M Q . 
This minimum mass is the same for both MFs (a) and (b) be- 
cause they are identical for to > 0.35M Q . Such a high mini- 
mal mass is not allowed in the disk, since many stars with mass 
to < 0.3 M Q are observed in the solar neighborhood. There- 
fore, either the bulge mass function is different from the disk 
MF, or the three long time scale events have a different origin 
(most likely stellar remnants). 

Without including the three afore-mentioned long-time 
events, which all correspond to a bulge clump giant source, the 
mean duration is t e = 16.7 days, yielding an average mass 
of 0.15 M Q . The minimum mass derived with the complete 
calculations to reproduce the time distribution in that case is 
n-inf = 0.11 Mq in case (a) and m- m { = 0.07 M Q in case (b). 
The KS test gives a result of 70% in both cases. Figure || shows 
the model distribution compared to observations in case (a). It 
is thus possible to explain the observed time distribution with a 
main sequence stellar population with the mass function (a) or 
(b) plus a stellar remnant population. 

4.1.2. Events whose source is a disk main sequence star 

The OGLE collaboration (Pacziriski et al., 1994) and the MA- 
CHO collaboration (Alcock et al., 1997), have observed the 
same region, near the Baade's window, in the direction of the 
bulge. The MACHO collaboration observed six events whose 
source lies on the disk MS, with a magnitude V ^ 18.5. The 
sources of 3 of these 6 events lie near the very red edge of the 
MS and could well be bulge blue helium core burning giants. 
This could be easily verified from their spectra. However, the 
remaining 3 sources are likely to be disk MS stars. These events 
yield an optical depth 7.5 x 10~ 7 , whereas the standard disk 
model (with no dark population) yields ~ 10~ 7 and a maximal 
disk (see paper II) yields ~ 2 x 10~ 7 . But the MACHO effi- 
ciency used for the derivation of this optical depth is an aver- 
age over all source stars. Since we are considering only sources 
with V < 18.5, the efficiency for these brighter sources must 
be higher. With an extreme (unlikely) 100% efficiency, the op- 
tical depth for disk sources is r = 2 x 10~ 7 in better agreement 
with disk models. We stress that theoretical estimates depend 
strongly on the disk scale length and scale height. 

The average time for the corresponding 3 events is (t e ) = 
4.4 days. Since the source stars are as bright as the bulge 
clump giants, the blending is expected to be negligible for these 
events. The mean duration would correspond to a minimum 
mass smaller than 0.005 M , definitely in the brown dwarf do- 
main. However, the extrapolation of either mass function (a) or 
(b) down to a minimum mass of 0.005 M requires a moderate 
power-law exponent in order not to exceed the local dynamical 
limit (Sbd = 27M pc~ 2 if a = 2 and E b d = 62M Q pc~ 2 
if a = 2.5 for MF (a)). Moreover, such a low minimum mass 
overproduces short time scale events for bulge sources, espe- 
cially when only clump giant sources are considered. 

Figure [| shows the calculated time distribution of the mi- 
crolensing events obtained with the MF (a) with two different 
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Fig. 4. Comparison of the giant bulge events (excluding the 3 
longest events) with three models. Dotted line: the MF (a) is 
assumed to be the same for both the bulge and the disk, with 
TOinf = 0.11 M Q . Dashed line: the disk MF has a m- m i = 
0.015 M and the bulge MF has a m inf = 0.4 M Q , for MF 
(a); the dot-dashed line is the same for MF (b). The histogram 
corresponds to observations. 

minimum masses for the bulge and the disk, as inferred from 
§4.1.1 and 4.1.2. A minimum mass m in f = 0.3 M Q for the 
bulge and m in f = 0.015 M Q for the disk, with both MFs (a) 
and (b). The model distributions agree with a 70% KS result 
with the observed one. Such a disk brown dwarf population 
corresponds to a local surface density: 

S bd w 20- 30 M pc- 2 (26) 

Clearly any solution between the ab ove o ne and the one 
where mj n f ~ 0.1 M everywhere (cf § |4.1.l| ), i.e. no brown 
dwarf, looks acceptable. The disk contribution to the bulge op- 
tical depth in this model is tj w 8 x 10~ 7 , slightly larger 
than the optical depth of a standard (no brown dwarf) disk 
(5 x 10~ 7 ), but still insufficient to explain the observed r w 
4xl0~ 6 . 

Whereas small number statistics and uncertainty in the dis- 
tance of the sources prevent a reliable determination of the 
brown dwarf contribution to the disk mass budget in this case, 
the present analysis corresponds to a brown dwarf and a dy- 
namical density in the disk, from (jl5|), (|T]) and (|26|): 

S M »0-30M oP r 2 ; E « 40- 7OM pc- 2 , (27) 

possibly larger than the value (fy determined in §2.1. A mass- 
dependent scale height does not change significantly the brown 
dwarf contribution, unless their scale height is larger than ~ 1 
kpc. 

These results leave open the following possibility: i) the 
mass functions of the disk and the bulge are different, in par- 



ticular in term of minimum mass, and ii) the dynamical up- 
per limit for S (Eq. ^ and the observed density of baryonic 
matter (Eq. [D]) imply that 20 to 30 M pc -2 of the solar sur- 
face density can be in the form of brown dwarfs, in agreement 
with the microlensing analysis (Eq. ^7|). Note that the possi- 
bility for the disk and bulge mass functions to be different 
has also been suggested by Gould et al. (1997) to reconcile 
star counts and microlensing observations. However, although 
these authors invoke a large population of brown dwarfs in the 
bulge, the present analysis suggest that they could as well be lo- 
cated in the disk. A more precise determination of the amount 
of brown dwarfs in the disk requires microlensing observations 
optimized for characteristic timescales of a few days, corre- 
sponding to a sampling of two observations per day. We stress 
the need for such microlensing searches. 

The present analysis illustrates the difficulty to reach robust 
conclusions about the mass under the form of sub-stellar ob- 
jects in the central regions of the Galaxy, and more precisely in 
the disk and the bulge, from present microlensing experiments. 
The situation will certainly be clarified once future projects will 
be operational, allowing much better statistics in both parts of 
the Galaxy. In fact, the MACHO survey has now a total of more 
than 200 candidate microlensing events, still under analysis, 
which will improve the present determinations. 

The volume density of disk brown dwarfs in the solar neigh- 
borhood deduced from the two types of calculations described 
above, global vs separate disk/bulge analysis, correspond re- 
spectively to pbd ~ 5 x 10 -3 M pc -3 (from Eq. p_l]) within 
about a factor 2, and pbd ~ 3 — 5 x 10~ 2 M Q pc~ 3 (from 
3 events in the disk, Eq. (|26|): this is clearly to be revised 
in the light of future observations), a factor 10 more. Future 
deep photometry and large field observational surveys might 
thus help determining the present issue by direct observation 
of nearby brown dwarfs, although the small density and the 
intrinsic faintness of brown dwarfs (My ^> 22, Baraffe et al. 
1997) renders the observation of a statistically-significant sam- 
ple of field brown dwarfs a tremendously difficult task. Note 
that infrared filters are highly recommended for these objects 
(Mj > 12; M K k, 11, Baraffe et al. 1997). 

The present analysis is thus an indication that an impor- 
tant amount of brown dwarfs in the disk is not excluded with 
present day observations, as suggested by the rise of the disk 
mass function near O.1M (cf. §3.1). 

4.2. Spheroid 

For a DeVaucouleurs spheroid, extrapolation of the MF 
into the brown dwarf region, even with m- m i = 0.01 M , yields 
a number of microlensing events towards the LMC of ~ 0.4 and 
an optical depth T sph — 5 x 10~ 9 (Chabrier & Mera, 1997). 
This illustrates the negligible contribution of the spheroid to 
the Galactic mass, and thus to the events observed towards the 
LMC. 
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4.3. Dark Halo 

The analysis of the first year of the EROS (Aubourg et al., 
1993) and MACHO (Alcock et al., 1993) microlensing obser- 
vations towards the LMC had shown that the observed events 
were likely to be due to halo brown dwarfs, with an average 
mass (m) « 0.03 M (Mera et al. 1996b; Kerins 1995). The 
inferred contribution of these objects to the missing mass was 
found to be between 10 and 20% (Alcock et al. 1995; Gates et 
al., 1995; Mera et al. 1996b). These calculations must now be 
re-examined in the context of the second year of the MACHO 
experiment. 

These results yield now a total of 1 to 2 events observed 
by EROS (Ansari et al., 1996; Renault et al., 1997) and 6 to 8 
by MACHO (Alcock et al., 1996). The inferred optical depth is 
much larger than derived previously, r b s = 2.2 ± 1 x 10~ 7 , 
about 40% of the value corresponding to the standard halo. The 
6 MACHO events^ yield (t e ) « 40 days, in reasonable agree- 
ment with the value inferred from the 2 events observed by 
EROS, (t e ) w 28 days. The statistics, however, are not suffi- 
cient to constrain a mass function from the time distribution. As 
will be shown in Paper II, the number of observed events can- 
not be explained by dark objects in the LMC itself. As shown 
by Chabrier & Mera (1997) and Graff & Freese (1996), extrap- 
olation of the MF ( |l6| ) into the brown dwarf domain yields a 
negligible fraction of brown dwarfs in the Galactic halo. 

A substantial amount of these objects in the halo thus im- 
plies a significantly different, much steeper MF beyond the 
hydrogen-burning limit. In order to estimate the maximum con- 
tribution of halo substellar objects to the microlensing counts, 
we consider, as in Mera et al (1996b), a halo MF ^(™) oc 
mT x with x = 5, essentially a Dirac -peaked MF. The maxi- 
mum normalization of §^ (to) at 0. 1 M Q is the observed count 
of halo subdwarfs in the solar neighborhood, derived in §3.4. 
We thus take as the upper limit for the dark halo substellar mass 
function: 



8 x 10 _3 M Q pc~ 3 is the local dark matter density (Bahcall 
1986), this dispersion corresponds to the asymptotic rotation 
velocity with the relation a = u ro t/v2 ~ 150 km.s -1 , where 
v TO t ~ 220 km.s -1 , as given by the isothermal sphere model 
(Binney and Tremaine 1987). This yields an average mass, cor- 
rected for blending, {\fm) 2 w 0.4 — 0.5 Mq, as suggested by 
the MACHO collaboration (Alcock et al., 1996). Several ob- 
servations (Beers & Sommer-Larsen, 1995; Dahn et al., 1995; 
Layden et al., 1996) suggest that the halo population velocity 
ellipsoid is radially elongated, by a factor ~ 1.5 — 2. This yields 
a projection effect which decreases the velocity dispersion. On 
the other hand, this dispersion is bound by the velocity of the 
line-of-sight toward the LMC (Eq. A. 19 ), which is (after proper 
calculations), (l/vios) -1 ~ 80 km.s~ . This corresponds to a 
minimum mass m m i = 0.04 M Q but to a zero dispersion ve- 
locity, clearly an irrealistic possibility. When considering such 
a "modified" halo model, with a non-isotropic velocity disper- 
sion tensor at ler of 80 km.s -1 and an important oblateness, 
the time distribution can be reproduced with mj n f ~ 0.1M Q , 
still above the hydrogen-burning limit. However, the number of 
observed events implies a normalization for the halo MF at 0. 1 
Mq significantly larger than the maximum value (|l6|), inferred 
from observed low-mass stars in the solar neighborhood. 

These calculations show that, if the results of the recent 
MACHO analysis are confirmed, the possibility for halo dark 
objects to be field brown dwarfs is clearly excluded, even for 
a MF substantially steeper than the ones in the disk and the 
spheroid, for any halo model. This analysis, consistent with 
the MACHO observations (Alcock et al., 1996) and with star 
count analysis at faint magnitudes (Chabrier & Mera, 1997), 
and the one conducted in §3.4 and 3.5, suggest a negligible 
mass-fraction under the form of field stars and brown dwarfs in 
the dark halo. Different solutions to try to reconcile star counts 
and microlensing observations in the halo will be considered in 
paper II. 



dJSl o 
— (to) = 4.0 x 10~ 3 
dm 



ni 



— 5 



,0.1M, 



l 0, 



(28) 

5. Summary and conclusions 



A Kolmogorov-Smirnov test done on the t e distribution of 
the EROS + MACHO halo events shows unambigously that 
such a MF is not consistent with the observations, whatever 
the value of m; n f. The optical depth can be reproduced only if 
minf < O.O2M , but the time distribution is in severe conflict 
with the observed one, as shown by a 0.02% KS probability 
((t e ) = 12 days). 

The expected average mass of microlensing events is re- 
lated to the characteristic timescale of the event and to the ve- 
locity dispersion (see Appendix Eqs. A2 and A15). The average 
tangential velocity is proportional to the dispersion of the Gaus- 
sian velocity distribution (l/v±)~ 1 oc a (see Appendix). As- 
suming a maximal dark halo ph(r) — Pdyn/r 2 , where payn = 



We have excluded the MACHO events #10, believed to be a vari- 
able star (Alcock et al., 1996), and #9, the binary event, which proba- 
bly belongs to the LMC. 



In this paper, we have examined in detail all the observational 
sources that constrain the mass distribution of the Galaxy. 
This includes dynamical constraints, circular rotation velocity 
and vertical velocity dispersion in the solar neighborhood, star 
counts in the disk, bulge and spheroid, and microlensing ob- 
servations towards the Galactic halo and the Galactic center. 
The star counts yield the derivation of stellar mass functions 
for the disk, the bulge and the spheroid down to the bottom 
of the main sequence. A substantial discrepancy, however, re- 
mains between the MF inferred from the nearby LF and the 
one deduced from the HST LF in the range 0.09 — 0.35 Mq. 
Although the nearby MF is consistent with a steadilly rising 
power-law dN/dm oc m~ a with a ~ 2 down to the bottom 
of the main sequence, the HST MF is essentially flat below 
0.6 M Q . Interestingly enough, the significantly rising last bin 
of the HST MF, indicates an upturn of the MF near the brown 
dwarf domain, so that both type of MFs, nearby and photo- 
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metric, suggest the presence of a substantial amount of brown 
dwarfs in the Galactic disk. 

We have determined the normalization of these mass func- 
tions and shown that the ratio of the halo to disk main sequence 
stellar density is about ~ 1/ 400 near the bottom of the main se- 
quence. The presence of the spheroid is proven unambigously 
by the observed star counts of high velocity stars but its con- 
tribution to the Galactic mass is negligible. These calculations 
yield the determination of the contribution of main sequence 
stars, and thus of observed baryonic matter, to the disk and halo 
mass budget. The detected (thin+thick) disk stellar surface den- 
sity corresponds to S„j s w 43 ± 5 M Q pc -2 , whereas the halo 
(spheroid + dark halo) density corresponds to less than ~ 1% 
of the dynamically determined dark matter, significantly less 
than the limit determined by the HST (Flynn et al., 1996). 

The contribution of sub-stellar objects to this mass depends 
on the lower limit down to which the mass functions can be 
extrapolated. This limit is fixed by the microlensing observa- 
tions, in particular by the time distribution of the events. This 
combined analysis of star counts for the luminous part of the 
mass function, which yields the determination of the slope and 
normalization near the brown dwarf limit, and of microlens- 
ing experiments, which determines the minimum mass, yields 
a consistent determination of the dark matter under the form of 
star-like objects, low-mass stars and brown dwarfs, in different 
parts of the Galaxy. Although uncertainties due to either small 
statistics or ill-constrained velocity dispersion and location of 
the sources in the bulge prevent for now a precise determina- 
tion, the derived brown dwarf surface density in the solar neigh- 
borhood, for a scale height h z = 320 pc (Gould et al., 1997), is 
likely to be in the range £bd = — 30 M Q .pc~ 2 so that the total 
local density in the disk is Sd yn ~ 40 — 70 M Q pc~ 2 . As will 
be shown in Paper II (Mera et al., 1997), this is in agreement 
with various observational determinations. More statistics of 
microlensing events towards the bulge will allow a separate de- 
termination of the brown dwarf contribution in the bulge and in 
the disk. An attempt based on the analysis of the events whose 
source is a clump giant or a disk main sequence star, respec- 
tively, yield rather surprising results (hampered by present ob- 
servational uncertainties): most brown dwarfs would be present 
in the disk, whereas the bulge would contain almost none of 
these objects. This stresses the need for further observations in 
this region. 

For the halo, we show that brown dwarfs have a negli- 
gible contribution (^ 1%) to the halo mass. A steep (a <; 
2) mass function in the halo seems to be excluded both by 
the microlensing analysis and by the M-dwarf LF. Therefore, 
the present analysis excludes main sequence stars and brown 
dwarfs as a significant contribution to the halo mass budget. 

The nature of the dark events observed in the halo remains a 
puzzle. Halo white dwarfs remain the least unlikely candidates, 
although this scenario implies severe constraints on the age of 
the halo and its initial mass function (Chabrier et al., 1996; 
Adams & Laughlin, 1996). More attention will be devoted to 
this problem in Paper II. 



These calculations, which combine observational con- 
straints arising from star counts, microlensing experiments and 
kinematic properties, yield the consistent determination of the 
amount of dark matter under the form of stellar and sub-stellar 
objects in the differents parts of the Galaxy. This yields new in- 
sight on the distribution of baryonic dark matter in the Galaxy 
and bears important consequences for the derivation of a con- 
sistent Galactic mass-model. This will be examined in the next 
paper (Mera, Chabrier and Schaeffer, 1997, paper II). 
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Appendix A: Microlensing equations 

The basic microlensing equations can be found for instance in Griest (1991) or in Kiraga and Paczyriski (1994). However, we use 
some extra formulae which are given below, after a summary of the basic equations. 

The duration t of a microlensingevent is defined as the time during which the amplification of the monitored star is larger than 
a given threshold amplification At\\ (usually At — 1.34). The event duration corresponds to the crossing time of the Einstein 
disk, of radius utR c , for the lens: 



1 = — v < - u — = 53 days v± ' x V^V 55k P c °- 5 v T " min (AA) 

where v± is the lens transverse velocity w.r.t. the line of sight, ut (resp. u m ; m ) is the impact parameter corresponding to the 
threshold amplification (resp. maximum amplification) At, L is the distance to the source, xL and m denote respectively the 
distance and the mass of the lens, and Re = -y/GmLx(l — x) is the Einstein radius. 
The characteristic time of an event is defined as: 



Re 



t e = ^ = —y/GLmx(l-x)= n/ . f , = (A.2) 



nun 



This effective time does not depend on the impact parameter, and can be recovered from the observations with the relation 
between u m i n and the known maximal amplification. In practice, the blending (several unresolved stars of which only one is 
amplified) renders this process rather difficult. The overall effect of blending is to underestimate the characteristic time. 

The probability for a source star to be microlensed at a given time is called the optical depth, and reads: 



u T 7ri? e L dx = / u T 7r — t. — p(xL)xCl — x) Ax (A3) 

m J c z 

where p(xL) = J °° p m (xL)dm is the total mass-density under the form of dark objects. 
The experimental optical depth is retrieved from observations by: 

_ Zt obs /e(t) 

Texp ~ N s xT obs {AA) 

where N s and T f, s denote the number of source stars and the total duration of the experiment, respectively whereas e(i) is the 
detection efficiency. 

• event rate: 

The theoretical event rate for a given Galactic model, i.e. the expected number of events per unit time, reads: 

dT = 2u T R e v ± ^^-P(m)P(v ± ) dm dv ± dx (A.5) 

m 

where P(m) and P(wj_) are the probability dis tribu tions respectivelly of lens mass and velocity. If the velocity distribution is 
independent of the position, the integration of (A.5) yields: 

4 fin™, i [L r + oo 

T = u T -VGL —=Y(m)dmx I y/x(l - x)p{xL)d{xL) I dw ± w ± P(w ± ) (A.6) 

c J rrlint v TO JO JO 

For N s source stars monitored during T b s , the expected number of events is: 
N = T x N s T obs (A.7) 



The relation (A.5) shows that the effective microlensing probability distributions are different from the model distributions, 
and are given by: 



eff 



(x) oc y/x(l - x)p(xL) (A.8) 



m 2 + 2 

7 or equivalentelly, a threshold impact parameter ut in units of the Einstein radius, related to At by At = 



UTy/u^ + 4 
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P eff(v±) OC V ± P(v±) 

P e// (m) oc 



in 

For a Maxwellian velocity distribution, of dispersion a, the transverse velocity reads: 
1 



¥{v) dv 
and 

Peff(v) 
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(A.9) 
(A. 10) 

(A.ll) 
(A.12) 



which yields: 

/ l^SV^ 1 - x)p{xL) dm dx 
and (with use of eqn. ( |A.6| )): 

32G 2 L 3 u T 



(A.13) 
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x / rnP(m) \x(l — x)\ p(xL)e c2cr2 'e dx dm 



(A. 14) 



For a given Galactic model, we use a Monte-Carlo algorithm to calculate the event rate, which yields a set of simulated events 
whose time distribution is exactly the one given by ( A. 14 ). 
These equations yield the average characteristic time (t e ): 



2VGL 



(A.15) 



In equation (A.15), the means are computed with the effective probability distributions (A. 8- A. 10). When writing explicitely 
the corresponding integrals, we can identify the optical depth and the event rate, which yields: 



t = r x -u T {t e ). 



(A. 16) 



A given Galactic model implies a time-distribution dT/dt e = Tx P(t e ) and the number of events predicted by the theory is 
given by: 



. . dr 

e(t e )—dt e 
dt e 



N th =Ex 

where E — N s x T b s 



(A.17) 



• Divergence of the moments: 

As shown in Eq. A.13, the probability to observe an effective time t e behaves as : 

P(t e )dt e oc \jt e P(v±)dv± oc v±dv± 



(A.18) 



so that P(t e ) has a long power-law tail at large times. This yields the divergence of the moments of order n > 3 J t™P(t e )dt e , 
which causes the extreme sensitivity of (i e ) and r to the rare large-time events. This long-range tail effect, which prevents the 
use of the higher order moments of t e , is also felt when the central limit theorem is invoked to identify the observed average time 
of e.g. the 45 events with the statistical average. They are equal in the large-N limit, but the convergence is extremely slow in 
this case. This shortcoming is fixed by using the average (j-) instead of (t e ). This method is more useful to constrain m ln i and 
(t e ) than Eq . A.15 since l/t e oc v± and is thus nearly gaussian at the beginning and l/t e oc m -1 / 2 and is thus most sensitive to 
the low mass end of the mass function. 
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• Motion of the line of sight: 

The lens velocity w.r.t. the line-of-sight (l.o.s.) depends on the motion of the l.o.s. itself and must include the motion of the 
Sun and of the source star w.r.t. the Galactic center. The lens velocity thus reads: 



V = VUnt. 



XV, 



(1 - x)Vq 



(A.19) 



The velocity now depends on the lens distance xL, and equations ( A. 1 1 ■ A. 14 ) relative to a gaussian velocity distribution, are 
no longer valid. Straightforward calculations yield the projected velocity of the lens w.r.t. the l.o.s.: 



(Vlent.,1 ~ V S< 1 + x(v Q j - V s j)) + (vi ent .,b ~ V Q . b + x(v Q . b - U a ,(,))' 



1/2 



(A.20) 



• Distance of the source star: 

In the case of the bulge, the elongation along the line of sight is no longer negligible, and the afore-mentioned formulae must 
include the possible variation of the distance L of the source, which implies an extra integral on L: 

= Iq +00 L^ s (L) dL^u 2 T ^p(xL)^xJT^x)dx 
f + °° L*v s (L) dL 

where v s is the density of source stars visible at the distance L. This requires a luminosity function for the source stars. Following 
Kiraga and Paczynski (1994), v s oc pL 2/3 , if the number of stars brighter than some absolute luminosity L is proportional to 
~Lr '. We have slightly modified v s to take into account the fact that giant stars have more or less the same luminosity. Then 
v s oc p(L 2/3 + C) where C is adjusted to reproduce the observed ratio of giants. 



The event rate (A. 6) now reads 



/ dm / dL / dx / d 2 vi ent / d 2 v s x -\d- v± v s (L)L 3 - 5 p(xL)y/x(l - x) (A.22) 

JO JO JO J-oo J-oa V m 



f +oo v s (L)L 2 dL 



This integral has six non-independent variables, only the mass can be separated. The Monte-Carlo integration method (see Press 
et al. 1992) does not rely on any discretization of the integration domain, and is hence suitable for this kind of integral. Moreover, 
the Monte-Carlo process provides a set of simulated microlensing events, from which the time distribution can be recovered 
easily. The adjunction of the experimental efficiency is also straightforward with a rejection algorithm. 



